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Abstract 

We study volume transition phenomenon in hydrogels within the framework of Flory-Rehner thermodynamic modelling; we show 
that starting from different models for the Flory parameter different conclusions can be achieved, in terms of admissible coexisting 
equilibria of the system. In particular, with explicit reference to a one-dimensional problem we establish the ranges of both 
temperature and traction which allow for the coexistence of a swollen and a shrunk phase. Through consideration of an augmented 
Flory-Rehner free-energy, which also accounts for the gradient of volume changes, we determine the position of the interface 
between the coexisting phases, and capture the connection profile between them. 
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1. Introduction 


Since the '80, temperature-induced discontinuous phase 
transitions were largely investigated from an experimental 
point of view in prolacrylamide gels and in nonionic N- 
isopropylacrilamide (NIPA) gels, either under free conditions 
or uniaxial tension ( [Tanaka et al.[ |1978| |1980| |1984[ |Hirotsu| 


[1987 Hirotsu and Onuki 1989| Suzuki and Kojima 1994 


[Suzuki et al. \991) . In particular, in (Hirotsu 1987| l, the results 
of the experiments were analyzed on the basis of the Flory- 
Rehner thermodynamic theory, in which the interactions be¬ 
tween polymer and solvent are accounted for by the Flory pa¬ 
rameter or polymer-solvent interaction parameter which mea¬ 
sures the dis-affinity between the polymer and the solvent. 
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When the Flory parameter is large, solvent molecules are ex¬ 
pelled from the gel and the gel shrinks, whereas when it is 
small, the gel swells (Doi 2009 j l. In spite of the large amount 
of data on temperature-dependent volume transition in gels, 
multiphysics-based models of thermally-driven volume tran¬ 


sition in hydrogels have only been recently proposed (Ji et al. 


[2006 1 [Birgesson et al.| [2008 [ [Duda et ah] [2010[ [Chester and[ 


Anand[[201 1[ Cai and Suo 201 1[ Ding et al. 2013[ Hong and 
Wangl [2013 Drozdov 2014 Drozdov and Sommer-Farsen 


[MT5] |. 

We make some progress in this topic by using an approach 
based on the Flory-Rehner free energy, appropriately extended 
to penalize gradients of volumetric strain. In this way we 
are able to describe coexistence between swollen and shrunk 
phases and to identify the interface position between them. In 
( Hong and Wang[ 2013| l, a similar approach was pursued, even 
if lacking of any predictions of the spatial distribution of the 
swollen and the shrunk phases at equilibrium. Our further con¬ 
tribution consists in discussing and solving the problem with 
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explicit reference to two different temperature-dependent mod¬ 
els of the Flory parameter: the first proposed in ([Hirotsu and 


Onuki 19891 for NIPA hydrogels and the second proposed in 
( |Afroze et aH |2000[ ) for aqueous solutions of uncrosslinked 
PNIPAM. It is worth noting that in (Hong and Wang |2013| l, 
only this last situation is considered, even if, as we explicitly 
show through a one-dimensional example, the two models de¬ 
liver different conclusions. 

With reference to a one-dimensional example, we propose an 
analysis, aimed to identify specific ranges of both temperature 
and traction which allow coexistence of the swollen and shrunk 
phases. The position of the interface is evaluated through an 
analytical method borrowed from (Cirillo et al. 2012| l, and 
shortly reviewed in the paper, under the limit for the higher or¬ 
der stiffness parameter, accounting for interfacial energy, going 
to zero. We found that the interface position between the co¬ 
existing phases depends on the temperature so that the relative 
portion of the sample occupied by the shrunk and the swollen 
phase change with temperature as well. We also verified this 
prediction through a numerical calculation. 


2. Equilibrium theory of gels 


we assume that both the volume of the solid and liquid compo¬ 
nent of the hydrogel do not change separately, and the change 
in volume of the system is a consequence of the variation of the 
fluid mass content, only; under this assumption, we have that 
J - 1/0. With this, the total energy density of the system only 
depends on the field J, and is parameterised by the temperature 
T of the hydrogel and the pressure P acting on the solid part 
of the system. Following ( |Doi 2009| l, and introducing the rel¬ 
ative change in volume between the dry and the current state, 
i.e., S = 7 - 1 > 0, we assume that the above mentioned total 
energy density G{S) per unit dry volume can be represented as 


RT^ 

r 7 




-H S log 


xS 


5 H- 1 5 -(- 




( 2 . 1 ) 


Here, we have: the gas constant R - 8.314 JK 'mol ', the 
solvent molar volume fl, the number of segments in a poly¬ 
meric chain, the scalar d accounting for the dimension of the en¬ 
vironment under consideration, and the Flory parameterbeing 
a temperature-dependent dimensionless parameter (called also 
the chi parameter), which represents the dis-affinity between 
the polymer and the solvent. 


Hydrogels are made of long-chain polymers which are 
cross-linked into a three-dimensional network and permeated 
by a solvent. From the point of view of continuum mechan¬ 
ics, they can be viewed as soft elastic materials consisting of an 
elastic matrix swollen with a fluid, and can be thermodynam¬ 
ically characterized through the choice of a free-energy. We 
start considering only isotropic deformation processes, which 
typically occur under free-swelling conditions; in this case, the 
state of the system can be described by the polymer fraction 0, 
which measures, locally, the fraction of volume occupied by the 
polymer, and by the change in volume J, which measures the 
change in volume from the dry state to the current state of the 
soft elastic material corresponding to the hydrogeQ Moreover, 


Let us note that the first addend in the equation (2.1 1 can be 
recognised as the isotropic elastic energy. Indeed, for d - 3 
(three-dimensional context) once introduced the shear modu¬ 
lus G = kbT Vc of the polymer, with kj, the Boltzmann constant 
and Vc the number of partial polymeric chains per unit volume, 
we get 3G((7'^^)^ - l)/2, which is exactly the isotropic elas¬ 
tic energy of three-dimensional elasticity when only isotropic 


deformation 7'^^ are involved (see ( Lucantonio et al. 2013 Q. 
For 7 = 2 and 7=1, analogous considerations can be made to 
identify the first addend with the isotropic elastic contribution: 
for 7 = 2, it holds under the assumption of plane deformations 
7'^^; for 7=1, under the assumption of uniaxial deformation 
7. 


*If we assume the material to be hyperelastic, the dependence of the free 


lated to the space dimension. The corresponding force is therefore given as the 


energy only on J is totally general, for isotropic deformations and is not re¬ 


derivative of the energy with respect to J. 
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2.1. Modeling the Flory parameterX 

One of the major assumptions of the Flory-Rehner thermo¬ 
dynamics model is that there is no volume change of the two 
components on mixing and that solvent and polymer can fit on 
the sites of the same lattice. It leads to a temperature indepen¬ 
dent additive constant in the expression of the Flory interaction 
parameter. Actually, these effects are not fully understood and 
all deviations from the lattice model, at the basis of the mixing 
component of Flory-Rehner free energy, are lumped into the in¬ 
teraction parameter x which can have non-trivial dependences 
on polymer fraction and temperature ( RubinsteTn} 2003[ ). Em¬ 
pirically, the temperature dependence of the Flory interaction 
parameter is written as the sum of two terms; the first is a tem¬ 
perature independent constant referred as the entropic part ofx', 
the second depends on temperature and is called the enthalpic 
part. With the aim to highlight the differences in phase tran¬ 
sition due to the model of x, we assume that the interaction 
parameter depends on the volume fraction 0 linearly, according 
to the model proposed in (Hirotsu and Onuki |1989| l for NIPA 
gels and in ( |Afroze et aH |2000| l for aqueous solutions of un- 
crosslinked PNIPAM. This choice can be thought of as the first 
order Taylor expansion of an unknown function of (ft: 


X^Xoj +X\.T^- 


( 2 . 2 ) 


As far as the dependence on the temperature is concerned we 
consider two different models. The first one, denoted from now 
on model a, prescribes that 


Xo.t-Aq + BqT and xu-^i+^iT, (2.3) 
with Ao, B(), A\, Bi reals; this is the choice which has been 


demonstrated experimentally to be valid for PNIPAM (Cai and 
Su"^ 120111) with Ao = 12.947, Bq = 0.04496K-', Aj = 17.92, 


and Bi - 0.0569K ^1. The second model, denoted from now 
on model b, prescribes 


Xo,t-Ao+^ and xij - M , 


(2.4) 


with Ao = 2.68294, Ai = 0.305, and Bq = 589.348K; this 
is the model proposed in (Hirotsu and Onuki 1989) 1 for NIPA 
hydrogels. 


3. Phase coexistence in hydrogels 

In this section we set the problem of phase coexistence for 
a hydrogel described by the energy density •HD- We start as¬ 
suming that a ,homogeneous phase of the system be a global 
minimum of the energy density ( |2.1| l corresponding to a con¬ 
stant field S and ask for the existence of parameters T and P of 
the model such that the system is in phase coexistence regime, 
that is to say, the energy density has multiple isolated local 
minima. 

We split the analysis into two cases: for the interaction pa¬ 
rameter X depending on the sole temperature T and for the in¬ 
teraction parameter depending on both the temperature T and 
the volume fraction 0; and discuss the nature of the solutions of 
the equation d Gjd S - 0. We show that forx depending on the 
sole temperature, the existence of multiple coexisting phases is 
possible only in dimension one {d = 1); whereas for;^f depend¬ 
ing on T and 0, it will turn out that coexistence is possible at 
any dimensions. 


3.1. Interaction parameter depending on the sole temperature 
We assume that the interaction parameter does not depend 
on the volume fraction. We then write x - Xt and the equation 
dGjdS = 0 for stationary states as 


RT. S 


{S + 1)5 


^-1 


1 


XT 


= -P. (3.5) 


’5-Hi 5-Hi (5 -h1)2J 

From now on the left hand side in p.5| ) will be denoted by L{S). 
In order to understand if equation ( |3.5| l admits one or more so¬ 
lutions, few remarks are needed on the behavior of L(S ) for the 
order parameter S tending to the admissible limits 0 and oo as 
well as on its minima. In particular we note that 


limL(5) = -co, (3.6) 

s-»o 

whereas 

lim L(S) = oo, (1/AJ(R7’/Q), 0, (3.7) 

S —»oo 

for d - 1,2,3, respectively. Moreover, we have that the equa¬ 
tion dL(S)ldS = 0 for the stationary points of L reads 

- 1 + (2xt - 1)5 = - 1 ) 5 (5 + 1)2/^+'. (3.8) 
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A graphical study of the above equation yields the following; 
for d -3 the function L has a single stationary point; for d - 2 
the function L has a single stationary point for xt >1/2 and no 
stationary point otherwise; for li = 1 there exists a real number 
1^ such that the function L has two stationary points for xt > 
XiNx), one stationary point for xt - Xi^x), and no stationary 
point otherwise (see figure [2 left panel). 


US) 





Figure 1: Left panel: qualitative study of equation (3.8^ in the case = 1; 
dotted and solid lines are the graphs of the functions on the left and right hand 
sides of the equation. Central & right panels: qualitative graphs of the function 
L(S ) introduced below 0’ XT'Xi^x) (center) and^r > x(Nx) (right). Note 
that for^ 7 '=;^(A^j:) an horizontal inflection point is present. 


one-dimensional case. This possibility is ruled out at larger 
dimensions. 


3.2. Free swelling state 


In the setup of the above Section 3.1 namely, under the as¬ 
sumption that the interaction parameter ^ 7 - only depends on T, 
we quickly recall the properties of the equilibrium state under 
free swelling conditions, i.e, P - Q (see also the discussion in 
(Doi |1996| l [Section 3.4]). From the result in the previous sec¬ 
tion, it is immediate that phase coexistence is not possible in the 
cases d - 2,3. We prove, here, that even ford= 1 the phase is 
unique. For d = I equation p.5[l becomes 


(5 + l)"log- 


1 1 


-—(S + iy-(S + l)^XTx (3.9) 


We denote by g(S ) the function at the left hand side of the above 
equation and note that g(S) -t-oo for 5 —> 0, g(S) -00 for 
S ^ CO, and 


From these remarks, it follows that the possible qualitative 
behavior of the function L are those depicted in figure [^(center 
and right panels). The structure of the equation ( |3.5] l suggests 
that it can have three solutions only for d - \ and xt > xi^x), 
and the energy function can have a double well structure only 
in such a case. Since for = 1 the energy G(S) tends to 0 and 
to -Hoo for 5 —> 0 and S ^ 00 , respectively, it follows, from 
the standard Maxwell equal area construction, that it is possible 
to find a value of the pressure/traction P such that the energy is 
characterized by two equally deep wells. 

In conclusion, in the case studied in this section, namely, for 
the interaction parameter not depending on the volume frac¬ 
tion (j), the hydrogel described by the total energy •HZ) can ex¬ 
hibit two coexisting homogeneous swelling states only in the 


^ It is possible to give a nice estimate of the number^(A^t). Indeed, for d = \ 
the equation 13.8) can be rewritten as 

—s'* + —5^ + —S^ + [ - (2xt - I) + —Is + 1 = 0 . 

Nx Nx Nx I ' Nx' 

Since, the signs of the coeflicient of the above polynomial can exhibit at most 
two variations, the number of its positive roots is at most equal to two. The 
condition for having two variations is> 1/2 + l/(2A^t). Hence we have that 
X(Nx)>i/2+l/(2Nxy 


^g(S) = 2(5 +1) log ^ ^ -1 - + If ■ (3.10) 

oS S S Nx 

Since the sum of the first two terms is smaller than or equal to 1, 

we have that dgjdS < 0 for any S. In conclusion, g is a function 

decreasing monotonically from +oo to -oo in the interval for 

S 6 [0, H-oo); hence, the equation p.9] l has a unique not trivial 

solution whatever the value of xt is. 

In the free swelling case {P - 0), the energy ( |2.1| l is such that 

G{S) 0for5 ^ 0andG(5) ^ +cx)for5 ^ +oo. Hence, the 

unique stationary point of G(5) is a minimum. In conclusion in 

the free swelling case, for any value of xtj namely, for any 

temperature, the system has a unique not trivial homogeneous 

phase. This result is in agreement with the discussion in ( |Doi[ 

[1996} . 

3.3. Interaction parameter depending on both the volume frac¬ 
tion and the temperature 


In this case, we get 

G(5) = 


RT) 

r d 


yWx 


^ 5 

+ S log + [xa.T 


(3.11) 
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Our original problem can be now rephrased as follows: are 
there any values of the thermodynamic parameters T and P 
which determine two different swelling states of the system? 
In other terms, we look for values of T and P such that the en¬ 
ergy G{S) has a double well graph with equally deep wells. For 
d = 2,3, it is possible to get a double well energy. However, 
we focus on the case d - being dimension one more suitable 
for setting and solving the interface location problem, through 
a technique already proposed in different contexts by some of 
the authors dCirillo et al.| [20091 |20T0l[20TTl[20T2l|20T3] l. 

Firstly, we note that S > 0 and 


lim G(S ) = 0 and 
s->o 


lim G(S) - +oo^ 

S—^OO 


Then, since the number of stationary points of the total energy is 
determined by its first partial derivative computed with respect 
to 5, we evaluate 




RT 1 1 


+Xl,T 


l-S 
{S + 1)3 


+ log 


S + 1 S + 


(3.12) 
— 1 -HP 

-H iJ 


and 


lim -—G(S) - —oo 
s->o dS 


and 


lim — G(5) = -Hoo^ 
5-^00 0 ^ 


The equation determining the stationary points dG{S)fdS - 0 
can be rewritten as 

T(5) = -P, 


(3.13) 


where 

L{S) 


RT 1 1 


Xu: 


l-S 


-Hlog- 


^1 


(3.14) 


(5 -H 1)3 ' 5 -H 1 S 

defines the traction exerted on the hydrogel. In order to estab¬ 
lish the number of solutions of the equation ( |3.13| l, we look at 
the graph of the function L(S). To do it, we note that 


lim L(S) - -oo and lim L(S) - -Hoo 

S^O s^oo 


(3.15) 


and 


±US) - 

dS Q. S{S + 1)4 La?;, 


La?, a?, a?. 


[—+ 1-2x0,T+2x1 j)s^ (3.16) 

{-ju + 2- 2x0,T - 4-xi,r)s + l] 


The numerator of dLjdS is a fifth order polynomial and its 
number of positive zeros can be estimated by means of the 
Descartes rule: first note that four of the six coefficients are 
positive and the two remaining are the first and the second or¬ 
der ones. We then have that, assuming that all the coefficients 
of the polynomial are different from zero, their sign can exhibit 
either zero or two variations. This implies that dLjdS has either 
zero or two positive real zeros. 

Hence, recalling ( |3.15| l, the graph of the function L{S) can 
be either monotonic or kinky (as we have already illustrated 
in figure center and right panel, in a different case)). The 
first case corresponds to absence of coexisting phases, whereas 
the second to the existence of two different phases. It is worth 
noting that the presence of two coexisting phase is possible only 
if the pressure is chosen properly via the Maxwell equal area 
rule. 

Now the natural question is: for what values of;^'o, 7 ' and;^^l J■ 
is the graph kinky? This question is not easy to answer, since 
studying the possibility of a quintic polynomial to have real 
roots is not trivial; we propose a nice estimate. With simple 
algebra we can write the derivative of the function L{S) as 


dS 


L(S) = 


-[SH1-2xo,t 


RT, I 1 

s(s + 1 ) 4 ' 

+2'Xi,t) (3-17) 

+ 25(1 -xo,T - 2x\,t) + l]}j 


and note that, a sufficient condition for dLldS to be strictly pos¬ 
itive is that the discriminant of the second order polynomial ap¬ 
pearing in the above equation is negative. Indeed, in such a case 
it would be also ensured that the first coefficient 1 -2xo,t+2x\,t 
is positive, so that the second degree polynomial is positive de¬ 
fined. Hence, we find the sufficient condition 


(1 -Xo,T - 2x\,t)^ - (1 - 2^0,r + 2xi,t) < 0, (3.18) 

which has the nice geometrical interpretation on the plane ^o,r- 
Xi,T depicted in figure]^ (light gray region). 

Another condition ensuring the derivative dLjdS to be pos¬ 
itive for 5 >0 is that the discriminant is positive but both the 
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Figure 2: The light gray region is a graphical representation of the solutions 
of the inequality i |3.18) . The dark gray region is a graphical representation 
of the solutions of the inequality i |3.19) . Note that the results depicted in this 
figure are consistent with those discussed in Section ED indeed, for ^ 1 ^ 7 - = 
0 the interaction parameter does not depend on the volume fraction, and the 
graph above is consistent with the fact that the system admits coexisting phases 
provided;^o,r is sufficiently large. 

coefficient of 5^ and S are positive too. Indeed, in such a case, 
by the Descartes’rule, it follows immediately that the two roots 
of the second order polynomial are both negative. Thus, we 
have the inequalities 

(1 -Xoj - - (1 - 2 ^ 0 , 7 ’ + ^-xu) > 0 

' ^ -Xa,T-2 xi,t > 0 (3.19) 

1 “ 2 ^ 0 ,r + 2.x\,t > 0 

which has the nice geometrical interpretation on the plane;^o, 7 — 
Xij depicted in figure[^(dark gray region). 

We remark that the conditions we found are sufficient, but 
not necessary, to rule out the existence of coexisting phases. 
Indeed, even if the second order polynomial in ( |3.17| l were neg¬ 
ative for some values of 5, it could happen that, due to the 
1/Vx additive term, the derivative dLjdS is positive for any S. 
Hence, the region with absence of multiple coexisting phases 
could be larger than the one depicted in figure We also note 
that such a condition becomes sharper and sharper when Nx is 
chosen larger and larger. 

4. Temperature driven phase transition in PNIPAM 

In this section we study the possibility of a temperature 
driven phase transition in PNIPAM in the framework of the 


one-dimensional model discussed in Section 13.31 We start as¬ 
suming it holds model a for the Flory interaction parameter. 
Under this assumption condition p.l8| ) reads 

(Bo + 2B,fT^ + 2[(Ao + 2Ai)(Bo + 2Bi) 

-3Bi]7’-h(AoH-2Ai)2-6Ai <0 (4.20) 

whereas conditions p.l9| l can be rephrased as follows 

(Bo + 2 Bi) 27’2 + 2[(Ao + 2Ai)(Bo + 2Bi) 

- 3Bi]r -H (Ao -H 2 Ai)2 - 6Ai > 0 

(4.21) 

1 - (Ao -H 2Ai) - (Bo -H 2Bi)7’ > 0 

1 - 2(Ao -H Ai) H- 2(Bo -H Bi)T > 0 

In the PNIPAM case ( |Cai and Suo||201 fj l, we have that ( |4.20| l 
is satisfied for 291.198 K< T < 301.868 K while ( |4~2T] ) is never 
satisfied. This means that the system can exhibit (depending on 
the pressure) two coexisting phases for T < Tiow < 291.198 K 
and T > Thigh > 301.868 K. We recall that the numerical esti¬ 
mate for Blow and Thigh are sharper and sharper provided Nx is 
larger and larger. 

This behavior is illustrated in in the first three panels of fig- 
ure[3 The Maxwell construction is illustrated at T = 283 K and 
r = 311 K in order to find the value of the coexisting pressure, 
namely, the pressure such that the two gray areas in figure 
or analogously the two local minima of the energy, are equal. 
Conversely at T = 300 K the energy has a single minimum, for 
any value of the pressure, so that a unique phase is observed. 

If the model b of the Flory interaction parameter is assumed, 
condition p.l8[ ) yields 

[-1 -h2Ao-2Ai -h(-1 +AQ + 2Ayf]T^ 

(4.22) 

h-2Bo(Ao-h2Ai)Th-B2 <0 
whereas conditions p.l9| l read 

[-1 -h2Ao-2Ai h-(-1 -hAo-h2Ai)2]t2 

-H 2Bo(Ao H-2Ai)T-H B^ > 0 

(4.23) 

-BoH-(l-Ao-2Ai)r>0 

-2Bo H- (1 -2Ao h-2Ai)T > 0 

We have that is satisfied for 126.869 K< T < 303.825 K 
while p.23| l is satisfied for T < 126.869 K. This means that the 
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Figure 3: Graphs of the traction L(S ) j3.14^ . First three panels with Flory parameter temperature dependence as in model a and, from the left to the right, T = 283, 
300, 311 K and P = 33.79, 15.00, -2.93 MPa. Fourth panel with FI 017 parameter temperature dependence as in model b,T = 315 K and P = -1.22 MPa. 


system can exhibit (depending on the pressure) two coexisting 
phases for T > Thigh > 303.825 K. This behavior is illustrated 
in figure (fourth panel). In the graph the Maxwell construc¬ 
tion is illustrated at T = 315 K and the coexistence pressure is 
found. 

We remark that the estimate for the Thigh coexistence temper¬ 
ature is very similar in the two theories. The main difference 
is in the fact that with the model a phase coexistence is possi¬ 
ble even at low temperature (T < Tiow), while with the model b 
coexistence is possible only at high temperatures. 



Figure 4: Coexistence pressures vs. temperature. Top: results for the model a 
at r < Flow are plotted. Bottom: results for the model a (black) and for the 
model b (gray) at T > Thigh are plotted. 


In the following we discuss the physical properties of the 
system in presence of two coexisting phases, the shrunk and 
the swollen one, which are the admissible equilibria with lower 


and higher volume change. 



320 340 360 380 400 

T[K\ 


Figure 5: Coexisting phases vs. temperature. Top: results for the model a at 
T < Flow are plotted. Bottom: results for the model a (black) and for the model 
b (gray) at T > Thigh are plotted. 


In figure]^ we plot the coexistence pressure versus the tem¬ 
perature. On the top we consider the model a at T < Tiow and 
note that the coexistence between the two phases occurs under 
compression. To the best of our knowledge the possibility of 
attaining coexistence in such regime has never been discussed 
in the related literature. On the bottom we consider both model 
a and model batT > Thigh, where we stress that the two models 
differ in the value of T^gh (see the discussion below ( |4.21| i and 
( |4.23| l). The coexistence between the two phases occurs here 
under tension for almost all the considered values of tempera¬ 
ture, but for the model a at temperatures close to Thigh, where 
coexistence occurs under compression. It is worth to notice that 
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in the T-P plane there exist a region, the one between the two 
curves, where the two models predict different phases for the 
system. 

In hgurej^we plot the coexisting phases versus the tempera¬ 
ture. On the top, we consider the model a at T < T\o^ and, on 
the bottom, both model a and model b alT > Thigh. It is worth 
to notice that the volume changes associated to the shrunk and 
the swollen phase are quite similar when T < T\o^, while they 
signihcantly differ one from the other at T > Thigh. Increasing 
the absolute difference between the current and the limit tem¬ 
perature, either Tiow or Thigh, the shrunk and the swollen phase 
tend to separate from each other. 


5. Interface location 


Gradient theories are suitable to be developed for model¬ 
ing stress/strain concentration due, for instance, to the pres¬ 
ence of geometrical singularities (crack propagation in fracture 
mechanics) or phase transitions as in the case of wetting. To 
the best of our knowledge, gradient theories have rarely been 
formulated to infer interface location when temperature-driven 
volume transition occurs in hydrogels. Some recent results 
within the framework of phase-held models have been obtained 
Hong and Wang ( 2013| l. 


in 


Within the one-dimensional framework introduced above, 
the gradient theory will be used as a tool to capture the posi¬ 
tion of the interface between two coexisting phases differing in 
the degree of swelling (shrunk and swollen phases). 

We consider the hnite interval [0,1] and introduce the Lan¬ 
dau energy functional 


r{S] 


-f 


G(S) + -(S'f 


dx 


(5.24) 


with the standard free-energy component G(S), dehned as in 
the equation ( |3.1 l| l, and a: > 0 an appropriate higher order 
stiffness; the prime denotes the derivative with respect to the 
space variable. We look for the equilibrium prohle S (x) by per¬ 
forming a standard variational computation and assuming either 
Dirichlet or homogeneous Neumann boundary conditions. We 


therefore get the following Euler-Lagrange equation 

dG 


kS" = 


dS' 


(5.25) 


We consider now Dirichlet boundary conditions corresponding 
to the shrunk (5sh) and the swollen (5sw) phases for the values 
of P and T ensuring coexistence. In other words, we solve the 
problem (|5.25|l with the boundary conditions 


5(0) = 5sh and S(l) = S, 


(5.26) 


By exploiting the one-dimensionality of the model a phase 
space analysis proves that the problem ( |5.25| l endowed with the 
above mentioned boundary conditions has a unique solution im¬ 
plicitly given by the integral 


r 


VA-di 


(5.27) 


Is., V2[£. + G(s)] 
where for any a- > 0 we have dehned implicitly by the equa¬ 
tion 

A/^d.s 

= 1 (5.28) 


I 


is., V2[£. + G(s)] 
namely, the integral in ( |5.27| l with S - S^v/ and x = 1. 

The prohle given by ( |5.27| l is a connection between the two 
phases and, in the limit a- —> 0 , presents a hat interface localized 
at 

V52G(5sw)/552 

X/ = , - - , ^ (5.29) 

V52G(5,h)/552 + ^Jd^GiS,JldS^ 

see (Theorem 2) in ( |Cirillo et al.||2012l l. We note that, if the en¬ 
ergy G had equal second derivatives at the phases, the interface 
would fall in the middle point of the interval. We underline that 
in our case the energy ( |2.1| l is not symmetric and therefore the 
interface position will change depending on the temperature T. 

In hgure we show that the interface tends to a sharp kink 
when the stiffness is chosen smaller and smaller. We introduce 
a physically more relevant parameter, the characteristic length 
{ dehned as 


f = - K. 

RT 


(5.30) 


We plot the solution of the equation ( |5.25| l with Dirichlet 
boundary conditions ( |5.261 l for different values of the charac¬ 
teristic length ( at given temperature and pressure ensuring co¬ 
existence. The plots show that the interface tends to localize in 
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Figure 6: Profile j5.29) for the model a with T = 283 K and P = 33.79 MPa (left), and T = 308 K and P = -1.78 MPa (center). Profile j5.29) for the model b with 
r = 308 K and P = -0.70 MPa (right). 



Figure 7: Interface position ^5.29^ as a function of temperature at coexisting 
pressure. On the left results for the model a at F < Flow are plotted. On the 
right results for the model a (black) and for the model b (gray) at F > Fhigh are 
plotted. 


the small stiffness limit. In figure (left and central panel) we 
solved the problem for the model a, whereas in figure (right 
panel) the model b was considered. We remark that the bound¬ 
ary value problem ( |5.25[ )-( |5.26| ) has been numerically solved 
by means of a properly implemented finite element code. 

It is worth noting that, as already remarked above, due to 
the fact that even at coexistence the energy G is not symmetric 
with respect to the central local maximum the position of the 
interface depends on the temperature. From the physical point 
of view, this means that at different temperatures the relative 
portions of the sample occupied by the shrunk and the swollen 
phases change. 

It is interesting to compare the center and right panels in fig¬ 
ure 1^ At temperature T = 308 K, the model a predicts that 
the interface position is close to the boundary where the shrunk 
phase is preserved; this means that the sample is mostly oc¬ 
cupied by the swollen phase. On the other hand, at the same 
temperature, the model b predictions are slightly different, in¬ 
deed the interface position is close to 0.3 so that a not negligible 
part of the system is occupied by the shrunk phase. 

In figure the interface position is plotted as a function of 
the temperature for both models. The difference between the 
predictions of model a and model b is pointed out in the picture 
on the right. From both pictures we conclude that at temper¬ 
atures far from the limiting values Tio* and Thigh the interface 
position is close to zero so that the sample is mostly occupied 
by the swollen phase. Conversely, the interface moves towards 
the center of the sample when the temperature gets closer to 
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its limiting values. This effect is much more important for the 
model b. 

We remark that our result goes in the opposite direction with 
respect to that illustrated in ( |Cai and Suo| |201 1| Fig. 14d). 
There, indeed, it is stated that at large temperature the shrunk 
phase tends to fill up the whole sample. However, we have to 


notice that the problems considered in our paper and in (Cai and 
Suo 2011| l are similar, but not equivalent, indeed, the setup in 
that paper is fully three-dimensional, whereas our discussion is 
limited to dimension one. 


6. Conclusions 

We addressed the modeling of thermally-driven volume tran¬ 
sition in hydrogels, within the Flory-Rehner thermodynamic 
setting; precisely, accounting for the temperature-depending 


pattern of the interaction parameter determined in (Hirotsu and 
Onuki 1989| l for NIPA hydrogels, and in (Afroze et al. 2000| l 
for aqueous solutions of uncrosslinked PNIPAM. In both mod¬ 
els, the Flory parameter depends linearly on the volume frac¬ 
tion of polymer within the gel; on the contrary, the dependence 
on temperature is different. We proposed a detailed analysis 
aimed to establish the ranges of both temperature and traction 
which allow for the coexistence of two different phases (swollen 
and shrunk) in the hydrogel. With specific reference to a one¬ 
dimensional problem, we showed as different models for the 
interaction parameter deliver different conclusions. 

Finally, for the values of temperature and traction ensuring 
phase coexistence, we presented a gradient model, appropri¬ 
ately extending the Flory-Rehner energy, and localize the in¬ 
terface position between the shrunk and the swollen phase, on 
one hand by means of a phase space analysis, previously devel¬ 
oped by some of the Authors, on the other developing proper 
finite element numerical calculations. 
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